Investigation of double beta decay with the NEMO-3 detector 
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Abstract 

The double beta decay experiment NEMO 3 has been taking data since February 2003. The aim 
of this experiment is to search for neutrinoless (0u/3/3) decay and investigate two neutrino double 
beta decay in seven different isotopically enriched samples ( 100 Mo, 82 Se, 48 Ca, 96 Zr, 116 Cd, 130 Te 
and 150 Nd). After analysis of the data corresponding to 3.75 y, no evidence for 0vf3f3 decay in the 
100 Mo and 82 Se samples was found. The halfdife limits at the 90% C.L. are 1.1 • 10 24 y and 3.6 • 10 23 
y, respectively. Additionally for 0vf3/3 decay the following limits at the 90% C.L. were obtained, 
> 1.3 • 10 22 y for 48 Ca, > 9.2 • 10 21 y for 96 Zr and > 1.8 • 10 22 y for 150 Nd. The 2v00 decay halfdife 
values were precisely measured for all investigated isotopes. 
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I. INTRODUCTION 



Interest in neutrinoless double-beta decay has seen a significant renewal in recent years 
after evidence for neutrino oscillations was obtained from the results of atmospheric, solar, 
reactor and accelerator neutrino experiments (see, for example, the discussions in [l-3]). 
These results are impressive proof that neutrinos have a non-zero mass. However, the 
experiments studying neutrino oscillations are not sensitive to the nature of the neutrino 
mass (Dirac or Majorana) and provide no information on the absolute scale of the neutrino 
masses, since such experiments are sensitive only to the difference of the masses, Am 2 . The 
detection and study of 0v(3(3 decay may clarify the following problems of neutrino physics 
(see discussions in {4 (J): (i) lepton number non-conservation, (ii) the nature of the neutrino 
(Dirac or Majorana particle), (iii) absolute neutrino mass scale (a measurement or a limit 
on rrii), (iv) the type of neutrino mass hierarchy (normal, inverted, or quasidegenerate), (v) 
CP violation in the lepton sector (measurement of the Majorana CP- violating phases). 

The currently running NEMO 3 experiment is devoted to the search for 0^/3/3 decay and to 
the accurate measurement of two neutrino double beta decay (2u/3/3 decay) by means of the 
direct detection of the two electrons. This tracking experiment, in contrast to experiments 
with 76 Ge, detects not only the total energy deposition, but other parameters of the process. 
These include the energy of the individual electrons, angle between them, and the coordinates 
of the event in the source plane. Since June of 2002, the NEMO 3 detector has operated 
in the Frejus Underground Laboratory (France) located at a depth of 4800 m w.e. Then in 
February 2003, after the final tuning of the experimental set-up, NEMO 3 has been taking 
data devoted to double beta decay studies. The first results with 100 Mo, 82 Se, 150 Nd and 
96 Zr were published in {7 12] . 



II. NEMO-3 EXPERIMENT 



The NEMO 3 detector has three main components, a foil consisting of different sources of 
double beta decay isotopes and copper, a tracker made of Geiger wire cells and a calorimeter 
made of scintillator blocks with PMT readout, surrounded by a solenoidal coil. The detector 
has the ability to discriminate between events of different types by positive identification of 
charged tracks and photons. A schematic view of the NEMO 3 detector is shown in Fig. 1. 
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NEMO 3 




The NEMO 3 detector is cylindrical in design and is composed of twenty equal sectors. 
The external dimensions of the detector with shields are 6 m in diameter and 4 m in height. 
NEMO 3 is_^ based on the techniques tested on previous incarnations of the experiment 
NEMO 1 



13|, and NEMO 2 Q. 



The wire chamber is made of 6180 open octagonal drift cells which operate in Geiger 
mode (Geiger cells). A gas mixture of ~ 95% helium, 4% ethyl-alcohol, 1% argon and 
0.15% water at 10 mbar above atmospheric pressure is used as the filling gas of the wire 
chamber. Each drift cell provides a three-dimensional measurement of the charged particle 
tracks by recording the drift time and the two plasma propagation times. The transverse 
position is determined from the drift time, while the longitudinal position is deduced from 
the difference between the plasma propagation times at both ends of the cathode wires. The 
average vertex position resolution for the two-electron events is at = 0.5 cm in the transverse 
plane of the detector and G\ = 0.8 cm in the longitudinal plane. The Geiger cell information 
is treated by the track reconstruction program based on the cellular automaton algorithm, 
described in 15]. 



The calorimeter, which surrounds the wire chamber, is composed of 1940 plastic scintil- 
lator blocks coupled by light-guides to very low-radioactivity photomultiplier tubes (PMTs) 
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TABLE I: Investigated isotopes with NEMO 3 [1 



Isotope 


100 Mo 


82 Se 


130 Te 


116 Cd 


150 Nd 


96 Zr 


48 Ca 


Enrichment, 


97 


97 


89 


93 


91 


57 


73 


% 
















Mass of 


6914 


932 


454 


405 


36.6 


9.4 


7.0 


isotope, g 

















developed by Hammamatsu. The energy resolution FWHM of the calorimeter ranges from 
14.1 to 17.6% for 1 MeV electrons, while the time resolution is 250 ps at 1 MeV. 

The apparatus accommodates almost 10 kg of different double beta decay isotopes (see 
Table 1). Most of these isotopes are highly enriched and all are shaped in the form of 
thin metallic or composite foils with a density of 30-60 mg/cm 2 . Three sectors are used for 
external background measurements and are equipped respectively with pure Cu (one sector, 
621 g) and natural Te (1.7 sectors, 614 g of nat TeO"2). Some of the sources, including 100 Mo, 
have been purified in order to reduce their content of 208 T1 (from 232 Th and 228 Th, and from 
228 Ra with a half-life of 5.75 y) and 214 Bi (from 226 Ra with a half-life of 1600 y) either by 

n □ 

a chemical procedure [16J, or by a physical procedure pj|. The foils were placed inside the 
wire chamber in the central vertical plane of each sector. The majority of the detector, 12 
sectors, notes a total of 6.9 kg of 100 Mo. 

The detector is surrounded by a solenoidal coil which generates a vertical magnetic field of 
25 Gauss inside the wire chamber. This magnetic field allows electron-positron identification 
by measuring the curvature of their tracks. The ambiguity of the e + /e~ recognition based 
on the curvature reconstruction is 3% at 1 MeV. 

The whole detector is covered by two types of shielding against external 7-rays and 
neutrons. The inner shield is made of 20 cm thick low-radioactivity iron which stops 7-rays 
and slow neutrons. The outer shield is comprised of tanks filled with borated water on the 
vertical walls and wood on the top and bottom designed to thermalize and capture neutrons. 

At the beginning of the experiment, the radon inside the tracking chamber, and more 
precisely its decay product, 214 Bi present in its radioactive chain, was found to be the 
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predominant background. Radon is present in the air of the laboratory and originates from 
the rock surrounding. It can penetrate the detector through small leaks. A tent coupled 
to a radon-free air factory was installed around the detector in October 2004 in order to 
decrease the presence of radon inside the tracker. 

Since February 2003, after the final tuning of the experimental set-up, NEMO 3 has 
routinely been taking data devoted to double beta decay studies. The calibration with 



radioactive sources is carried out every six wee 



daily with a laser based calibration system 17 1. 



cs. The stability of the calorimeter is checked 



The advantage of the NEMO 3 detector rests in its capability to identify the two electrons 
from [3/3 decay and the de-excitation photons from the excited state of the daughter nucleus. 
The NEMO 3 calorimeter also measures the detection time of the particles. The use of 
appropriate time-of-ffight (TOF) cuts, in addition to energy cuts, allows for an efficient 
reduction of all backgrounds. 

A full caption of the detector and its characteristics can oe found in Q. 

III. EXPERIMENTAL RESULTS 

The f3(3 events are selected by requiring two reconstructed electron tracks with a curvature 
corresponding to a negative charge and originating from a common vertex in the source foil. 
The energy of each electron measured in the calorimeter should be greater than 200 keV. 
Each track must hit a separate scintillator block. No trackless PMT signals are allowed. The 
event is recognized as internal using the measured time difference of the two PMT signals 
compared to the estimated time-of-flight difference of the electrons. 

The background can be classified into three groups: external 7-rays, radon inside the 
tracking volume and internal radioactive contamination of the source. All three were esti- 



mated from NEMO 3 data with events of various topologies [18[. In particular the radon 
and internal 214 Bi is measured with eja events. The ej, e'yy and 6777 events are used to 
measure the 208 T1 activity requiring the detection of the 2.65 MeV 7-ray typical of the 208 T1 
/3-decay. Single electron events are used to measure the foil contamination by /3-emitters. 
The external background is measured with the events with the detected incoming 7-rays 
that produce an electron in the source foil. The control of the external background is done 
with two-electron events originating from the pure Copper and natural Tellurium foils. Mea- 
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TABLE II: Two neutrino half-life values for different nuclei obtained in the NEMO 3 experiment 
(for 116 Cd, Ca and 130 Te the results are preliminary). The first error is statistical and the second 
is systematic while S/ B is the signal-to-background ratio. 



Isotope 


Measurement 
time, days 


Number of 
2v events 


S/B 


2i/ 2 (2z/), y 


i nrin r 

100 Mo 

100 Mo- 100 Ru(0+) 

82 Se 

116 Cd 
96 Zr 

iso Nd 

48 Ca 
iso Te 


389 
334.3 

389 
168.4 
1221 

939 
943.16 
1152 


219000 
37.5 
2750 
1371 
428 
2018 
116 
236 


40 
4 
4 

7.5 
1 

2.8 
6.8 
0.35 


(7.11 ± 0.02 ± 0.54) x 10 18 [8] 
(5.7±J:g ± 0.8) x 10 20 [10] 
(9.6 ± 0.3 ± 1.0) x 10 19 [8] 
(2.8 ±0.1 ±0.3) x 10 19 
(2.35 ± 0.14 ± 0.19) x 10 19 [12] 
(9.11±g;H ± 0.63) x 10 18 [11] 
(4.4+°;^ ± 0.4) x 10 19 
(6.9 ± 0.9 ± 1.0) x 10 20 



surements performed with HPGe detectors and with radon detectors are used to verify the 
results. Different level of radon activity in the tracking volume were detected in the data 
taken before October 2004 (Phase 1), and after the installation of the antiradon facility 
(Phase 2). 



A. Measurement of 2v(3(3 half-lives 

Measurements of the 2vf3(3 decay half-lives have been performed for seven isotopes avail- 
able in NEMO 3 (see Table 1). The NEMO 3 results of 2v/3j3 half-life measurements are given 
in Table 2. For all the isotopes the energy sum spectrum, single-electron energy spectrum 
and angular distribution were measured. The 100 Mo double beta decay to the Of excited 
state of 100 Ru T 2 / 2 = [5.7±H(stat) ± 0.8{syst)} ■ 10 20 y has also been measured by NEMO 



3 10]. For 100 Mo, 82 Se, 96 Zr and 150 Nd these results are published. For the other isotopes 



their status is preliminary. 
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B. Search for decay 



No evidence for Qu/3/3 decay was found for all seven isotopes. The associated limits are 
presented in Table 3. 

TABLE III: Limits at 90% C.L. on Oz^/3/3 decay (neutrino mass mechanism) for different nuclei 
obtained in the NEMO 3 experiment. 



Isotope Measurement Ti/ 2 (0^), y 



time, days 


100 Mo 


1409 


> 1.1 


•10 24 


y 


82 Se 


1409 


> 3.6 


• 10 23 


y 


i30 Te 


1221 


> 1 • 


10 23 y 


150 Nd 


939 


> 1.8 


.1 22 


y 


116 Cd 


77 


> 1.6 


• 10 22 


y 


48 Ca 


943.16 


> 1.3 


.1 22 


y 


96 Zr 


1221 


> 9.2 


•10 21 


y 



The Oi//3/5-decay search in NEMO 3 is most promising with 100 Mo and 82 Se because of the 
larger available sample mass and high enough Qpp ~ 3 MeV. The data taken from February 
2003 till the end of 2008 has been used in the search for neutrinoless double beta decay. 
The data corresponds to 1409 effective days of data collection. For the case of the mass 
mechanism, the Oi//3/3-decay signal is expected to be a peak in the energy sum distribution 
at the position of the transition energy Qpp. The two-electron energy sum spectra for 100 Mo 
and 82 Se are shown in Fig. 2 and Fig. 3 demonstrating a good agreement between observed 
and expected number of events. 

Since no excess is observed at the tail of the distributions, limits are set on the neutrinoless 
double beta decay with the CL S method based on the use of a log-likelihood ratio (LLR) 



test statistics 



19(]. This gives the lower limits on the half-life of T^ 2 ( 100 Mo) > 1.1 ■ 10 2 V 



(90% C.L.) and T^ 2 ( B2 Se) > 3.6- 10 2 V (90% C.L.) A lower half-life limit translates into an 
upper limit on the effective Majorana neutrino mass (m u ). According to the recent results 
of different theoretical NME calculations (see Table 4) the half-life limit obtained for 100 Mo 
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FIG. 2: Distribution of the energy sum of two electrons for 100 Mo (left) and 82 Se (right), 1409 d 
data. The shape of a hypothetical 0u signal is shown by the curve in arbitrary units. 
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FIG. 3: Distribution of the energy sum of two electrons in the region around Q/3/3 value for 100 Mo 
(left) and 82 Se (right), 1409 d data. High energy tail of the energy sum distribution for events 
in molybdenum (left) and selenium (right) foils are shown with black points. The background 
contributions are shown within the histogram. The shape of a hypothetical 0v signal is shown by 
the curve in arbitrary units. 



corresponds to the neutrino mass interval (m u ) < 0.45 - 0.93 eV. It is less restrictive for 
82 Se (m u ) < 0.89 - 2.43 eV. The reached NEMO 3 sensitivity on the neutrino mass is close 



to that of IGEX 



20|, the Heidelberg-Moscow collaboration 21] and CUORICINO [22J. 
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TABLE IV: Limits on the effective neutrino mass {m u ) (in eV) corresponding to different theoretical 
model calculations of nuclear matrix elements obtained for T]/ 2 (Ov/3/3) > 1.1 • 10 24 y in case of 
100 Mo and T 1/2 (0z//3/3) > 3.6 • 10 23 y in case of 82 Se. 



Nuclear matrix elements 100 Mo 82 Se 



Shell model [23] - < 2.43 

QRPA [24], [25] < 0.58 - 0.75 < 1.12 - 1.38 

QRPA [26] < 0.45 - 0.93 < 0.89 - 1.61 

IBM-2 [27] < 0.49 - 0.55 < 1.03 - 1.19 

PHFB [28] < 0.70 



C. Search for double beta decay with Majoron emission 

The 0h>x°(3l3 decay requires the existence of a Majoron. It is a massless Goldstone boson 
that arises due to a global breakdown of (B-L) symmetry, where B and L are, respectively, 
the baryon and the lepton number. The Majoron, if it exists, could play a significant 
role in the history of the early Universe and in the evolution of stars. A 2/3-decay model 
that involves the emission of two Majorons was proposed within supersymmetric theories 
and several other models of the Majoron were proposed in the 1990s (see review [29] and 
references therein). The possible two electrons energy spectra for different 0vx°l3(3 decay 
modes of 100 Mo are shown in Fig. 4. Here n is the spectral index, which defines the shape 
of the spectrum. For example, for an ordinary Majoron n — 1, for 2v decay n = 5, in the 
case of a bulk Majoron n = 2 and for the process with two Majoron emission n = 3 or 7. 

No evidence for Oz^x°/3/3 decay was found for all seven isotopes. The limits for 100 Mo, 
82 Se, 150 Nd and 96 Zr are presented in Table 5. In particular, strong limits on "ordinary" 
Majoron (spectral index 1) decay of 100 Mo (T 1/2 > 2.7 • 10 22 y) and 82 Se (T 1/2 > 1.5 • 10 22 
y) have been obtained. Corresponding bounds on the Majoron-neutrino coupling constant 
are < g ee >< (0.35 — 0.85) • 10~ 4 and < (0.6 — 1.9) • 10~ 4 , respectively (using nuclear matrix 



elements from 



23-123). 



9 



250 - 



200 



150 - 



100 



50 



2^XX 

/ \n=5 2Px n=1 
i=2 




500 1000 1500 2000 2500 3000 

Energy, KeV 



FIG. 4: Energy spectra of different modes of 2(32v (n = 5), 2/3%° (n = 1 , 2 and 3) and 2f3x°X°( n 
3 and 7) decays of 100 Mo (see text). 



TABLE V: NEMO 3 limits on Ty 2 (y) for decay with one and two Majorons at 90% C.L. for modes 
with spectral index n = 1, n = 2, n = 3 and n = 7. 



Isotope 


n =1 


n = 2 


n = 3 


n = 7 


100 Mo [8] 


>2.7-10 22 


>1.7-10 22 


>1-10 22 


>7-10 19 


82 Se [8] 


M.5-10 22 


>6-10 21 


>3.1-10 21 


>5-10 20 


96 Zr [12] 


>1.9-10 21 


>9.9-10 20 


>5.8-10 20 


>1.1-10 20 


150 Nd [11] 


M.5-10 21 


>5.4-10 20 


>2.2-10 20 


>4.7-10 19 



IV. CONCLUSION 

The NEMO 3 detector has been operating within its target performance specifications 
since February 2003. The 2v(3f3 decay has been measured for seven isotopes with high 
statistics and greater precision than previously. The 100 Mo 2v/3(3 decay to the 0^ excited 
state of 100 Ru has also been measured. No evidence for Oz^/3/3 decay was found for all seven 
isotopes. The best 0u/3/3 limits have been obtained with 100 Mo (> 1.1 • 10 24 y at 90% C.L.) 
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and 82 Se (> 3.6 • 10 23 y at 90% C.L.). 

NEMO 3 is current experiment with data collection continued into ~ 2011. 
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